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I. INTRODUCTION 

In (1), the single space scatterer formalism based on (2) is used to evaluate the bulk 
parameters of a cloud illuminated from outside. The single space scatterer formalism 
allows an obstacle excited by an incident wave traveling in free space to radiate in free 
space. The calculations of (1) took only into account contributions due to dry air and 
water. The approximations for the attenuation agree closely with those obtained using 
number density and the scattering cross section of water drops. The results are obtained 
directly and they represent a solid starting point for subsequent developments. 

In (3), we used a modified two-space scatterer formalism of Twersky (4) and (5) to 
establish for a cloud modeled as a statistically homogeneous distribution of spherical water 
droplets, the dispersion relations that determine its bulk propagation numbers Ky and bulk 
indexes of refraction ly ( j— 1,2) in terms of the vector equivalent scattering am plitude Q 
and the dyadic scattering amplitude g of the single water droplet in isolation. The results 
of (3) were specialized to the forward direction of scattering while demanding that the 
scatterers preserve in the sense of (5), the incident polarization. This requirement did 
allow us in (3) to drop the subscript j and to look for K and 77. 

Here, we apply (3) to obtain specific numerical values for the macroscopic parameters 
of the cloud. We work as in (6), with a cloud of density p— 100 cm“ 3 , a wavelength 
A = 0.7774 //m, and with spherical water droplets of common radius a= 10 pm. In addition, 
the scattering medium is divided into three parts, the medium outside the cloud with 
propagation parameter k Q = ^ — 8.0823 ^m” 1 , moist air ( the medium inside the 
cloud but outside the droplets) corresponding to k x - k^, and the medium inside 
the spherical water droplets specified by ac 2 = KiT} 2 - The physical constants rj 1 , and 
*?2 — 1 -33+ 2 7-33 x 10 -8 (from (7)) are the relative indexes of refraction of moist air and 
water respectively. The numerical values of K and rj are important to the innovative work 
of (8) for lightning due to point sources inside the cloud. 

In this modified two-space scatterer formalism, a single spherical water droplet is 
allowed to be excited by a wave traveling in K- space (the space of the equivalent medium) 
but to radiate in - space instead of the usual free or k 0 - space. This is possible since we 
have a sparse distribution of scatterers ( i.e., the particles in the distribution are located at 
a significant distance away from each other as compared to their radius) and the boundary 
correction terms due to back and forth reflections from k 1 — space to k 0 — space or vice 
versa are negligible. 

The results of this report are applicable to a cloud of any geometry since the boundary 
does not interfere with the calculations. Also, it is important to notice the plane wave 
nature of the incidence wave £ in the moist atmosphere. Should the cloud be of a 
geometry such that its direct response to an initial outside excitation be of spherical 
type, one would have to make the plane wave approximation in order to apply the work 
of Twersky. The theory of Twersky is mostly suitable for a slab distribution of scatterers. 
However, it can be extended to cover the present problem without a lost of generality since 
the contributions of the boundary layer surfaces are small and most lightning activities 
usually occur inside the cloud. 
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In order to avoid repetition, we cite key equations of (3—5, 10 ). In general, we work in 
spherical coordinates. We use bold face or an arrow to denote a vector or a vector operator. 
A circumflex indicates a vector of unit magnitude. A tilde on the top of a letter denotes 
a dyadic (second rank tensor). For brevity, we use [5:4] for equation 4 of Re. ( 5) etc. 

II. MATHEMATICAL ANALYSIS 

For an incident plane harmonic electromagnetic wave <f> x = a 1 e , (* l ‘ r_,w< ) in the moist 
air due to an initial outside excitation <j> 0 of the cloud, the self-consistent integral equation 
governing the multiple configurational scattering amplitude as in (5) is 

G <( f ) = g<(f, Ki) • a. 1 e' K I ‘ + J g t (r, r c ) ■ G m (r c )e' g “ Rtm , [1] 

C 

where Rtm = r< — r m , f ~ ^ f dfl c . The single dyadic scattering amplitude g(r, &*) • a t = 

c 

g(r, k x : a x ) is defined in (3:4). The magnitude of the separation distance |R im | < D ( the 
diameter of the cloud). 

After taking the ensemble average of [1], using the quasi- crystalline approximation of 
Lax (9), the equivalent medium approach, and Green’s theorems, the dispersion relations 
determining the bulk parameters (for more details see (5)) are 

G (** I = J~ V { [ e "'^' R - }+P f [/( R ) - [2] 

< H K >- K V Vb — w 

where dR denotes volume integration over (Vj> — v). Here, Q is the equivalent scattering 
amplitude and Uy is a radiative function defined by 0 j — J g(r , r c )-Q ^/c lc | e*' ?lc ' R , and 

c 0 = «i/47rz. The bulk propagation parameter is K j = k x rfj with iy being the bulk index of 

refraction, and {[Z,^]} = / [fd n g — gd n f]dS is the Green surface operator. The ensemble 

S 

is specified as in (5) by the average number p of scatterers in unit volume and by pf( R) 
with /(R) as the distribution function for separation of R pairs. 

Similar to (3), the model is required to neglect all phase transition effects (10) and to 
take only into account pair interaction due to central forces. Neglecting the inter-droplet 
potential, the distribution function /( R) can be chosen to be always equal to unity. Hence, 
[2] is simplified to 

(K 2 - <)t + (^)g(f, £)].<?(*, | K) = 0. [3] 


In [3], we let r — K and use the fact that the scatterers preserve the incident polar- 
ization to transform [3] into 


,K- - «;> . -(*£),(*,*), „■ - . - -(*?>(«,*) 


[4] 
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where the subscript j is no longer necessary. In [4], the two-space scatterer formalism 
scattering amplitude corresponding to the spherical droplet, a large tenuous scatterer, 
must be evaluated in the forward direction of scattering (i.e., g = g(r, r ) = a a . g^K, K^. 

The unit vector K, the direction of the bulk propagation vector K, is such that K • a 2 = 0. 


III. NUMERICAL CALCULATIONS 

To start the numerical calculations, we use the WKB approximation corresponding 
to a bulk excitation associated with an incident electric field tp = xe* K2r for the two- 
space scatterer formalism scattering amplitude g = g(r, f) = x-g(K, k) in the forward 
direction. The WKB approximation (12) consists of replacing the field inside the scatterer 
by a field traveling in the direction of the incident excitation which propagates inside the 
scatterer with propagation constant of the medium of the scatterer. For the two space 
scatterer formalism, the propagation parameter of the incident wave is different from that 
of the radiated wave. Therefore, the conventional WKB method is not applicable here. 
Using a modified version given in (4) and the two space scatterer formalism volume integral 
representation of the scattered wave, it can be shown that 


Here, as in (4) and (12), 


[5] 
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Inserting [5] into [4] yields the desired results for the bulk propagation n umb er K 
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and the bulk index of refraction rj 
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where the physical constants /i 2 , and p are the magnetic permeability of the medium 
inside the water droplet and the bulk magnetic permeablity respectively. The formulae for 
the calculation of //, and e in the forward direction of scattering are obtained from (11:28) 

= 1 - 2 Ttpi^F — T*) /tjkI, e = i + 2 irp(j r -f f*) //c 2 , 

IK X 


To work with [7-9], we need to calculate the value of k x = k 0 r\ x = k 0 (rj xr -f- tjy lf *), and 
k 2 = k x tj 2 — k x (1.33+27.33 x 10 The real part r) 1T of r) x is obtained from the formula 
given by (13) 


(?? ir - i)io 6 = [0.378,125 + 

p{i+(i-Q 49 ~ o.oi5t)p x io ~ 6 } 
1+0.003661* 


0.002,141,4 

0.000,017,931 

A 2 H 

A4 J 


0.0624 — 


0.000,680 


[ 10 ] 


A 2 J 1+0.003661* 


where from (14), f (the vapor pressure at *— io°C) = 9,209 mmHg , and the pressure 
p — 720 mmHg is the lowest permissible for [10]. With A = 0.7774/2772 in [10], we have 
7 ] xr — 1.00026476, and K xr = 8.08443. 

To obtain the imaginary part k xx of k x , we use the result of (15) for the attenuation 
of intensity in the air Ie~ ad — 0.8/, where I is the intensity, a = 2 K a i ( the subscript 
a is used to indicate air), and the path is taken to be 11 km. From the above, we have 
Kai = 1-014 x 10 -11 . Using Archimedes principle, the value of « a ;, and the sphericity of 
the water droplets in the cloud, we can write k x = 8.08443 + *1.0446 X 10” 11 , « 2 = 

10.75229 + z.59260 x io“ 6 . 

To proceed with these calculations, we must verify the validity of the approximation 
that gives [5]. The WKB approximation is valid according to (12) if the host medium 
is such that 3?(/c 1 )2a(e ir — 1) >> 1 and (e ir — 1) < 1 where e x = e ir + ie x j, e ir = 
Vir ~ Viii e ii — 2 VirVii since the magnetic permeability p x is assumed to be one. In this 
case , the above conditions are satisfied. From the imaginary part of tc x , we deduce the 
value of 7 j x i = ( K x i/k 0 ) = 1.2924 x io~ 12 . 

Equations [7], [8], and [9] form a coupled nonlinear system for the bulk parameters K, 
77, 6, and //. To solve the system numerically, we approximate the exponential factor of l x 
by its leading term and use the values given in (1) as initial data. Hence, we have 


K = 8.08442974208285 + *3.14117832893745 x 10“ 8 , 
e = 1.00000000014550 + *3.87231835597287 x 10 -9 , 
p = 1.00000000014570 +*3.89603100094467 x 10 -9 , 

77 = .999999968400000 + *3.88417455500000 x 10“ 9 , 
T = 15.1452985315156+ *404.033753935163, 

F* = -1.01021251084684 x 10 -2 + *1.23330047641423. 


The results of [11] are close as expected, to those given in (1). Should the incident wave 
were magnetic in nature, e and p would have been interchanged. To improve [11], one 
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can up grade the model and construct /( R) numerically. The values given in [11] can be 
directly applied to the problem of lightning due to point sources located inside the cloud. 
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